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Abstract. Using full general relativistic calculations, we investigate the possibility of generation of mass outflow 
from spherical accretion onto non-rotating black holes. Introducing a relativistic hadronic-pressure-supported 
steady, standing, spherically-symmetric shock surface around a Schwarzschild black hole as the effective physical 
barrier that may be responsible for the generation of spherical wind, we calculate the mass outflow rate Rm 
in terms of three accretion parameters and one outflow parameter by simultaneously solving the set of general 
relativistic hydrodynamic equations describing spherically symmetric, transonic, polytropic accretion and wind 
around a Schwarzschild black hole. Not only do we provide a sufficiently plausible estimation of Rm, we also 
successfully study the dependence and variation of this rate on various physical parameters governing the flow. 
Our calculation indicates that independent of initial boundary conditions, the baryonic matter content of this 
shock-generated wind always correlates with post-shock flow temperature. 
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^ 1 1. Introduction 

>t i Accretion onto massive celestial bodies is believed to be responsible for energy generation mechanism in many astro- 
physical objects including galactic binary systems and extra-galactic objects like Quasars and AGNs. Investigation 
of transonic, spherically symmetric accretion onto gravitating astrophysical objects was carried out from a purely 
Newtonian approach as well as using a pseudo-general relativistic and complete general relativistic framework (Hoyle 
& Littleton 1939, Bondi 1952, Michel 1972, Shapiro 1973a,b, Blumenthal & Mathews 1976, Begelman 1978, Brinkmann 
1980, Malec 1999, Das 2001, Das & Sarkar 2001, Sarkar & Das 2001). However, in works mentioned above, spherical 
accretion was shown to exhibit essentially directed motion (radial inflow) which smoothly passes through a sonic 
point and no discontinuity in accretion onto a black hole was taken into account. Nevertheless, it was understood 
that accretion onto a black hole could be efficient in producing radiative energy if the directed in-fall motion could 
be randomized close to the hole by some suitable mechanism the best of which would be shock formation by some 
means. For radial accretion, a novel mechanism has been proposed (Protheros & Kazanas 1983, Kazanas & Ellision 
1986, hereafter PK83 and KE86 respectively) in which the kinetic energy of spherically accreting material has been 
randomized by proposing a steady, collision-less, relativistic hadronic-pressure-supported spherically symmetric shock 
around a Schwarzschild black hole which produces a nonthermal spectrum of high energy (relativistic) protons. A 
thermal distribution of particles may also be produced at the shock which radiates by thermal bremstrahlung but only 
the radiation resulting from the non-thermal (which might actually contain most of the energy, see Axford 1981) was 
considered in their model. Since radial in-fall (freely falling situation) in a gravitational potential is unlikely to release 
amounts of energy sufficient to produce the observed high luminosity and high energy tail of the spectrum for most of 
the quasars, it has been suggested in PK83 and KE86 that the relativistic particles required to account for the observed 
radio, X-ray and 7-ray emission may result from shock acceleration of a part of in-falling plasma via first order Fermi 
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acceleration which transforms directed kinetic energy into relativistic particle energy and thus produces relativistic 
protons with high efficiency A steady-state situation is assumed to have developed where a standing collision-less 
spherically symmetric shock is proposed to form. 

The question of the existence of a steady shock for spherical accretion deserve some elaborate discussion and clarifi- 
cation. As the absence of considerable angular momentum in accreting material does not supports standard Rankinc- 
Hugoniot type shock formation, assumption of a steady and standing spherically-symmetric shock for freely falling 
material, may apparently contradict the concept of shock formation in general fluid flow. The situation, however, is dif- 
ferent if the relativistic protons of sufficient energy density are present around the black hole, ft is worth studying what 
happens when freely falling supersonic (supersonic, because, as we will see later, spherical accretion with a reasonably 
realistic value of inflow energy and accretion rate crosses its sonic point far away, ^ hundreds of Schwarzschild radius 
or even more, from the black hole) plasma approaches close (some factor of tens of Schwarzschild radius or even less) 
to the black hole. Due to various possible instabilities in supersonic plasma flow, any small disruption mechanism may 
slow down the flow and create a 'piston' which produces a shock. Thus the initial formation of shocks may be attributed 
either to convection of ambient relativistic particles or to dissipation in magnetic field as suggested in the literature 
(see PK83, KE86, and references therein). Once the shock is produced, it accelerates a part of the in-falling matter 
to relativistic energies. Infalling thermal particles with free fall velocity are then assumed to be shock-accelerated via 
first order Fermi acceleration (Axfrod, Leer & Skardon 1977, Blandford & Ostriker 1978, Blandford & Eichler 1987, 
and references therein, Kirk & Duffy, 1999, and references therein) and relativistic protons will be produced. These 
shock accelerated relativistic protons were assumed to be trapped within a radius comparable to the radius of the 
shock surface by some mechanism (possibly by small amount of stochastic magnetic fields tied to the accreting matter, 
which is small enough to make the hydrodynamical approximation valid to describe the flow) because had it been the 
case that those protons would not be confined to the central region of the quasars, spallation of nuclei in the outer gas 
cloud responsible for the line emission in Quasars would result in an overabundance of boron (Baldwin et. al. 1977), 
contrary to the observation. In addition, the matter density in this region may be assumed to be sufficiently high so 
that the relativistic protons are depleted mainly by nuclear interactions which, possibly could set a lower limit (though 
not very restrictive) on the accretion rate for a given black hole mass (PK83). It was also suggested in PK83 that 
if the particles accelerated at the shock were neucleons rather than electrons, then, due to their negligible Compton 
losses, the Compton catastrophe could be avoided. 

Because of the fact that the kinetic energy of these relativistic particles are much larger than their gravitational poten- 
tial energy, and because their scattering mean free path could be larger than the black hole radius, [] these relativistic 
protons suffer practically no Compton losses and are not instantly swallowed by the black hole, rather they scatter 
several times and undergo multiple inelastic nuclear collisions and produce pions (7r ±: °) before being captured by the 
black hole (PK83, KE86). Pions generated by this process decay into relativistic electrons, positrons, neutrinos and 
antincutrinos and produce high energy 7-rays (see §2 of Das (1999) and references therein for various particle produc- 
tion and decay mechanism) . As the relativistic protons are assumed to lose energy only through the above mentioned 
nuclear collisions and since the proton-proton collision time scale t pp would, in general, be much larger compared to 
the free-fall time scale, a sufficiently high density of relativistic protons could be maintained which provides radially 
outward pressure sufficient enough to support the sustenance of a steady and standing shock. The shock is, thus, 
self-supported. 

In a previous work (Das 1999), we have explicitly calculated the exact location (radial distance measured from the 
central accretor in units of Schwarzschild radius r g = 2GA C { BH , where Mbh is the mass of the black hole, G is the 
Universal gravitational constant and c is the velocity of light in vacuum) of the above mentioned steady, standing, 
spherically-symmetric shock in terms of only three accretion parameters, namely, the specific energy of the flow, accre- 
tion rate (scaled in units of Eddington rate) and the adiabatic index of the flow for spherically-symmetric, transonic 
accretion of adiabatic fluid onto a Schwarzschild black hole. By solving the set of hydrodynamic equations describing 
the motion of accreting material under the influence of pseudo-Schwarzschild potential proposed by Paczyhski & Wiita 
(1980), we have shown there that it is possible to construct a self-consistent inflow-outflow system where a part of the 
accreting material may be blown as spherical wind from the spherical shock surface and the mass outflow rate Rm (the 
measure of the fraction of accreting material being 'kicked off' as wind) was computed in terms of various accretion 
as well as shock parameters. However, calculation presented in our previous work was based on the valid assumption 
that spherically symmetric, transonic flow onto a Schwarzschild black hole can be well approximated using a suitable 
modified Newtonian potential and no general relativistic calculations were performed to study the flow profile. 
In this paper, we would like to construct a self-consistent, spherically symmetric, polytropic, transonic, non-magnetized 
accretion- wind system using a full general relativistic description of flow structure around a Schwarzschild black hole. 
In doing so, we would like to perform the following steps: 



Furthermore, the converging magnetic field, if present, is also expected to mirror the relativistic particles and help prevent 
them from being absorbed by the black hole, see KE86. 
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(a) We will write down the set of general relativistic hydrodynamic equations of motions describing the accretion and 
wind in Schwarzschild metric. 

(b) Incorporating a suitable equation of state, we will solve the equations to get a general flow profile for accreting 
material where any dynamical (flow velocity Mach number of the flow etc.) as well as thermodynamic (flow temper- 
ature, density, pressure etc.) quantity can be explicitly calculated at any radial distance (measured from the accretor 
in units of r g ) in terms of only three accretion parameters, namely, the total specific energy and adiabatic index of 
the flow and accretion rate scaled in units of Eddington rate. 

(c) We will use the formula for shock location obtained in Das (1999) with the assumption that direct use of that 
formula in this paper would not be quite unjustified because procedure to calculate the shock location in modified 
Newtonian description would roughly be unchanged also in the full general relativistic framework. This is because the 
shock location is obtained by equating the two alternative expressions for luminosities of the same object, hence any 
correction term arising from the red-shift (gravitational or cosmological) would cancel out. Using the values of dynam- 
ical (and thermodynamic) variables computed in general relativistic framework in terms of three accretion parameters 
mentioned in (b), we will calculate the shock standoff distance. As the condition necessary for the development and 
sustenance of such a shock is satisfied when at shock location, the Mach number of the inflow M s h is considerably 
high (high supersonic flow at the time it encounters the shock, Ellison & Eichler 1984), we concentrate on polytropic 
accretion with such a set of values of conserved specific energy £ and accretion rate (scaled in the unit of Eddington 
rate) MEdd which produces a high shock Mach number solution. 

(d) We will take this hadronic pressure-supported shock surface to be the generating surface of mass outflow. At the 
shock surface, density of the post-shock material as well as the post-shock flow temperature shoots up and dynamical 
flow velocity falls down abruptly. In other words, highly supersonic inflow becomes subsonic and accreting matter 
becomes shock-compressed at this surface. Matter starts getting hotter and denser and starts piling up on the shock 
surface. The post-shock relativistic hadronic pressure and thermal pressure (pressure generated by high temperature 
produced at the shock) then gives a kick to the pilcd-up matter, the result of which is the ejection of outflow from 
the shock surface. Thus a hot and dense spherical shock surface serves as the 'effective' physical atmosphere regarding 
the generation of mass outflow from matter accreting onto black holes. At the shock, entropy is generated and the 
outflow as well as post-shock inflow will have higher entropy density for the same specific energy. For this type of 
inflow, accretion is expected to proceed smoothly after a shock transition, since successful subsonic solutions have 
been constructed for accretion onto black holes embedded within normal stars with the boundary condition u = c 
(Flammang 1982); where u is the in- fall velocity of matter and c is the velocity of light in vacuum. 

It is to be mentioned here that one fundamental criterion for formation of hydrodynamic outflows is that the out- 
flowing wind should have a positive Bernoulli constant which means that the matter in the post-shock region is able 
to escape to infinity. However, positiveness in Bernoulli's constant may lead to another situation as well where shock 
may quasi-periodically originate at some certain radius and propagate outwards without formation of outflows. So 
formation of outflows is one of the possible scenarios when we focus on the positive energy solutions. In this paper 
we concentrate only on solutions producing outflows, thus we use only the positive Bernoulli constant throughout our 
model which is standard practice for studying hydrodynamic winds. Another assumption made in this paper is to 
treat the accreting as well as post-shock matter as a single temperature fluid, the temperature of which is basically 
characterized by proton temperature. 

The plan of this paper is as follows. In next section we describe how to formulate and solve the governing equations. 
In §3, we present our results. In §4, the results will be reviewed and conclusions will be drawn. 

2. Governing equations and solution procedure 

2.1. Inflow model 

We assume that a Schwarzschild type black hole spherically accretes fluid obeying a polytropic equation of state. The 
density of the fluid is p(r), r being the radial distance measured in the unit of Schwarzschild radius r g . We also assume 
that the accretion rate (in the unit of Eddington rate MEdd) is not a function of r and we ignore the self-gravity of the 
flow. For simplicity of calculation, we choose the geometric unit where the unit of length is scaled in units of r gi units 
of velocity in units of c . All other physically relevant quantities can be expressed likewise. We also set G = c = 1 in 
the system of units used here. 
For a Schwarzschild metric of the form 




the energy momentum tensor T^ v for a perfect fluid can be written as 
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where e and p are proper energy density and pressure of the fluid (evaluated in the local inertial rest frame of the 
fluid) respectively and w M is the four velocity commonly known as 

= — 

ds 

Equations of motion which are to be solved for our purpose are, 
1) Conservation of mass flux or baryon number conservation: 



(puf,). M = 



(3a) 



and 

2) Conservation of momentum or energy flux (general relativistic Euler equation obtained by taking the four divergence 
of TV): 

(e + p) u KV u v = -p^ - Upp^u" (36) 
where the semicolons denote the covariant derivatives. 

Following Michel (1972), one can rewrite Eq. 3(a) and Eq. 3(b) for spherically symmetric flow as 

inpur 2 = M m 



and 



p + e 



1 



(4a) 



(46) 



as two fundamental conservation equations for time-independent hydrodynamical flow of matter on to a Schwarzschild 
black hole without back-reaction of the flow on to the metric itself. M,„ being the mass accretion rate and C is some 
constant (related to the total enthalpy influx) to be evaluated for a specific equation of state. 
For polytropic equation of state i.e, 

p = Kp lin 



(where ji n is the adiabatic index of the inflow) and defining M. (where M. = M in ^ in n*- - 1 K - 1 is a measure of 
the specific entropy of the flow) to be another conserved quantity of the flow except at the shock location (because 
shock generates entropy), one can rewrite conservation equations 4(a) and 4(b) as (Das & Sarkar 2001, and references 
therein): 



hu t 



1>- 



1 - ± 

r_ 

l-v? 



and 



M = 47T 



1) 



uutr 



(5a) 



(56) 



L7m - (a 2 + 1)_ 

where £ , h,u t and a are the total conserved specific energy( which include the rest mass of matter), specific enthalpy 
(h = T Tl(^2+i) ) i specific binding energy and local adiabatic sound speed respectively. Eqn. 5(b) may be considered 
as the outcome of the conservation of mass and entropy along the flow line. The expression for the adiabatic sound 
speed a can be written as (Weinberg 1972, Frank et. al. 1992): 



dp\ 2 



linP 



Hi 



,.kT 



(6) 



where T is the flow temperature, /x is the mean molecular weight and rriH~m p is the mass of the hydrogen atom and 
k is Boltzmann's constant. The subscript S indicates that differentiation is performed at constant specific entropy. 
One can now easily derive the expression for velocity gradient (^) (by differentiating Eq. 5(a) and 5(b)) as 



dfu _ u (1 - u 2 ) [a 2 (4r - 3) - l] 
dr ~ 2r(r- 1) (u 2 - a 2 ) 



(7a) 



Since the flow is assumed to be smooth everywhere, if the denominator of Eq. 7(a) vanishes at any radial distance r, 
the numerator must also vanish there to maintain the continuity of the flow. One therefore arrives at the so called 
'sonic point (alternately, the 'critical point') conditions' by simultaneously making numerator and denominator of Eq. 
7(a) equal to zero and the sonic point conditions can be expressed as follows 
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For a specific value of £ and ji n , location of sonic point r c can be obtained by solving the following equation algebraically 

64r c 3 (£ 2 - 1) + 16r c 2 (2£ \\ n - 9) + 4r c {£\\ n - 27) + 27 = (7c) 

where y in =(%to±f). 

To determine the behaviour of the solution near the sonic point, one needs to evaluate the value of (^) at that point 
(we denote it by (^) ) by applying L 'Hospitals' rule on Eq. 7(a). It is easy to show that (^p) can be obtained by 
solving the following quadratic equations algebraically: 

du\ 2 (7m - 1) (lgrg - 16r c - 8 7m r c + 6 7m + 3) /du\ 
rfr / c 3r c (4r c -3)^ V dr / C 

( 7m -l)(2r c -l)(24r c 2 -28r c -8r2 7m +4r c7m + 37 m + 6) 

+ 1 = (7d) 

2r c (4r c -3) 2 (r c - 1) 

It is now quite straightforward to simultaneously solve Eq. 5(a) and Eq. 5(b) to get the integral curves of the flow 
(curves showing variation of Mach number of the flow with radial distance) for a fixed value of £ and j in . Detail 
methodology for this purpose will be discussed in §2.3. 

In this work we normally use the value of 7 i„ to be |. Though far away from the black hole, optically thin accreting 
plasma may not be treated as ultra-relativistic fluid (by the term 'ultra-relativistic' and 'purely non-relativistic' we 
mean a flow with 7 = | and 7 = | respectively, according to the terminology used in Frank et. al. 1992) in general, close 
to the hole it always advects with enormously large radial velocity and could be well-approximated as ultra-relativistic 
flow. As because our main region of interest, the shock formation zone, always lies close to the black hole (a few tens 
of r g away from the hole or sometimes even less, see results and figures in §3), we believe that it is fairly justifiable to 
assign the value | for 7 j„ in our work. However, to rigorously model a real flow without any assumption, a variable 
polytropic index having proper functional dependence on radial distance (i.e., 7 = 7(r)) might be considered instead 
of using a constant ji n , and equations of motion might be formulated accordingly, which we did not attempt here for 
the sake of simplicity. Nevertheless, we keep our option open for values of 7j„ other than | as well and investigated the 
outflow for an range of values of "f in for a specific value of £ and M^dd (Fig- 5a, §3.3.3). The same kind of investigations 
could be performed for a variety of values of £ and Mem and a set of results may be obtained with a wide range of 
values of 7i„ which tells that our calculation is not restricted to the value ji n — | only; rather the model is general 
enough to deal with all possible value of % n for polytropic accretion. 



2.2. Shock parameters and the outflow model 

As already mentioned, a steady, collision-less shock forms due to the instabilities in the plasma flow. We assume that 
for our model, the effective thickness of the shock A s h is small enough compared to the shock standoff distance r s h , and 
that the relativistic particles encounter a full shock compression ratio while crossing the shock. At the shock, density 
of matter will shoot up and inflow velocity will drop abruptly. If (p_,u_) and (p + ,u+) are the pre- and post-shock 
densities and velocities respectively at the shock surface, then 

£±-B 

— ^comp — \° ) 

P- u + 

where R comp is the shock compression ratio. For high shock Mach number solution, the expression for R CO m P can be 
well approximated as 

Rcomp = lAAM sh i (9) 

where M s h is the shock Mach number and Eq. (5) holds for M s h £ 4.0 (Ellison & Eichler 1985). 
In terms of various accretion parameters, shock location can be computed as (Das 1999): 

_ SoppMsdd ( I - 2AM sh -°- 6S \ 

where <j vv is the collision cross section for relativistic protons, u s h and M s h are the dynamical flow velocity and the 
Mach number attained at the shock location, MEdd is the mass accretion rate scaled in units of Eddington rate. One 
can understand that 

(u sh ,M sh ) = C (£,MEdd,lin) (10a) 
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where £ has some complicated non-linear functional form which cannot be evaluated analytically, but the value of u a h 
and M s h can easily be obtained in terms of j<£, MEdd,7in} by numerically solving Eq. 5(a-b) and Eq. 10 (with the 
help of Eq. 7(a-d)) simultaneously. Hence one can write 



r s h 



ti(£,M E dd,lin) (106) 



where £ has some functional form other than that of £. 

In ordinary stellar mass-loss computations (Tarafder 1988, and references therein), the outflow is assumed to be 
isothermal till the sonic point. This assumption is probably justified, since copious photons from the stellar atmosphere 
deposit momenta on the slowly outgoing and expanding outflow and possibly make the flow close to isothermal. This 
need not be the case for outflow from black hole candidates. Our effective boundary layer, being close to the black 
hole, are very hot and most of the photons emitted may be swallowed by the black hole itself instead of coming out 
of the region and depositing momentum onto the outflow. Thus, the outflow could be cooler than the isothermal flow 
in our case. We choose polytropic outflow with a different polytropic index 7 Q < 7j„ due to momentum deposition. 
Nevertheless, it may be advisable to study the isothermal outflow to find the behaviour of the extreme case. Modelling 
the isothermal outflow is in progress and will be presented elsewhere. 

In our calculation we assume that essentially the post-shock fluid pressure and the post-shock proton temperature 
controls the wind formation as well as the barionic matter content of the outflow. The presence of a collision-less 
steady standing spherical shock discussed in this work may randomize the directed in-fall motion and at the shock 
surface the individual components of the total energy of the flow (which is a combination of the kinetic, thermal and 
gravitational energy) gets rearranged in such a manner that the thermal energy of the post-shock matter dominates 
(due to enormous shock generated post-shock proton temperature) over the gravitational attraction of the accretor 
and a part of the in-falling material is driven by thermal pressure to escape to infinity as wind. In §3.3 we show that for 
any shock solution, whatever the initial boundary condition obeyed by the pre-shock flow, the mass-loss rate always 
co-relates with the post-shock proton temperature which essentially supports the validity of our assumption. 
The adiabatic post-shock sound speed d^ h and the post-shock temperature T^ h (which is basically the temperature of 
the protons according to our one-temperature fluid approximation) can be calculated as: 



V Psh 



and 



T+ = ^ (11.) 



K Psh 

where p^ h and p^ h are the post-shock pressure and density of the flow at shock location r s h respectively. For low 
energy accretion ('cold' inflow, so to say) which is appropriate to produce a high shock Mach number solution, one 
can assume that the pre-shock thermal pressure {p~ h ) may be neglected compared to its post-shock value (pj/J and to 

the pre-shock ram pressure (yp~ h (u~ h ) 2 ^J . One can obtain the value of p^ h using Eq. (8-9) and from the total pressure 



balance condition at shock as, 



Pt h =W h ) 2 rs h (^P^) (12) 



Rcomp 



Combining Eq. (8-12), post-shock sound velocity and temperature obtained at the shock surface can be rewritten as: 

a f h = o.eoeiiw^ V 5 KJ- 1 - 5 (1 - omAM sh -^) - 5 Q = ^ (is«) 



and, 

^m p M E dd<T 3 p Ju+ h )~ 3 \ nent „ - 2.4M S/ T°' 68N " 



T+ = 0.014 (^ M ^y( u s h ) (1 _ . 694Ms?r o. 75) (i_ 



3.2M sh -° S2 



(136) 



The two conservation equations governing the outflow will exactly be the same as Eq. (5a) and (5b) with different 
polytropic index, i.e., with 7 = 7o for outflow. So we can write, 



£ = hu t (14a) 

and 

(7o-lK 



M = 47T 



7o - {al + 1) 



to) 2 

u u t r z (146) 
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where £ is the specific energy (including rest mass) of the outflow which is also assumed to be constant throughout 
• ' • i' i 

the flow and M = M out j ' ,0 ~ 1 K° is the entropy accretion rate of the outflow. 7o < 7i „ as already mentioned. 
Also Eq. (4a) can now be written as the mass conservation equation for the outflow as: 



M out = Att P °u ( r y 



(15) 



In Eq. (15), p° and u° are the density and dynamical flow velocity of the outflowing material. We then define the mass 
outflow rate R m as: 

M out 



Rri 



(16) 



Like Eq. (7a-7d), the value of sonic point conditions and the value of ^f- at the outflow sonic point can be written 



as: 



du° 
dr 



u° 


i - K) 2 ~ 




(a°) 2 (4r - 3) 


- 1 


2r (r - 1) 


V) 2 - (a°f 





1 



4r° - 3 



(4r c °) 3 



(V) 2 -l| +(4r c °) 2 [2 (£')\l-o\ +4r c ° \(s')\ 2 -27 



+ 27 = 



(17a) 

(176) 
(17c) 



where cp = {^fj. 



dr) c 3r°(4r°-3) f V dr 

( 7o - 1) (2r c ° - 1) [24 {r°f - 28r° c - 8 (r c °) 2 7o + 4r° 7o + 3 7o + 6 



+ - 



= 



2r°{Ar° c -2,)* (r° - 1) 
Where the superscript o indicates that the quantities are measured for the outflows. 



(17d) 



2.3. Simultaneous solution of inflow-outflow equations 

In this work, we are interested in finding the ratio of M out to M in (Eq. (16)), and not the explicit value of M out . Also 
note that the primary goal of our present work was to compute the outflow rate and to investigate its dependence on 
various inflow parameters but not to study the collimation procedure of the outflow. 

Before we proceed in detail, a general understanding of the transonic inflow outflow system in the present case is 
essential to understand the basic scheme of the solution procedure. Let us consider the transonic accretion first. 
Infalling matter becomes supersonic after crossing a saddle- type sonic point, the location of which is determined by 
j<£, Msdd, Tin j- This supersonic flow then encounters a shock (if present), location of which (r s h) is determined from 
Eq. (10). At the shock surface, part of the incoming matter, having a higher entropy density (because shock in a fluid 
flow generates entropy), is likely to return as wind through a sonic point other than the point through which it just 
entered. Thus a combination of transonic topologies, one for the inflow and other for the outflow (passing through a 
different sonic point and following topology completely different that of the 'self- wind' of the accretion), is required to 
obtain a full solution. So it turns out that finding a complete set of self-consistent inflow outflow solutions ultimately 
boils down to locate the sonic point of the polytropic outflow and the mass flux through it. Thus a supply of parameters 
£,MEdd, lin and 7o make a self-consistent computation of Rm possible. Here 7o is supplied as free parameter because 
the self-consistent computation of 7o directly using £, MEdd and 7in has not been attempted in this work; instead we 
put a constrain that 7o < 7 j„ always and for any value of 7in . In reality, 7o is directly related to the heating and 
cooling processes taking place in the outflow. 

We obtain the inflow sonic point r c by solving Eq. (7c). Using the fourth order the Runge Kutta method, u(r), a(r) 



and the inflow Mach number 



u(r) 
a(r) 



are computed along the inflow from the inflow sonic point r c till the position 

where the shock forms. The shock location is calculated by solving Eq. (10). Various shock parameters (i.e., density, 
pressure etc at the shock surface) are then computed self-consistently. 

For outflow, with the known value of £' and 7o , it is easy to compute the location of the outflow sonic point r° from 
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Fig. 1: Variation of shock location r sh (in units of r g ) with conserved specific energy £ (which includes rest mass energy) and 
accretion rate MEdd (scaled in units of Eddington rate) of the inflow for a fixed value of inflow adiabatic index ji n (=|). While 
£ and MEdd are plotted along X and Y axes respectively, r s h is plotted along the Z axis. It is observed that r s n correlates with 
M E dd and anti-correlates with £. 

Eq. (17c). At the outflow sonic point, the outflow velocity u° and polytropic sound velocity a° is computed from Eq. 
(17b). Using Eq. (17a) and (17d), and {^£) c is computed as was done for the inflow. Runge -Kutta method is 

then employed to integrate from the outflow sonic point r° towards the black hole to find out the outflow velocity u° 
and density p° at the shock location. The mass outflow rate is then computed using Eq.(16). 
It is obvious from the above discussion that Rm should have some complicated non-linear functional dependence on 
the following accretion and shock parameters: 

Rm = * (£,MEdd,r s h,M sh ,R comp ,^ in ,^ (18a) 
As r s h, M s h and R comp can be found in terms of £, MEdd and j in only, ultimately it turns out that: 

R m = ft (£, M E dd, 7m , 7o ) (186) 
Where f2 has some complicated functional form which cannot be evaluated analytically. 

3. Results 

3.1. Shock Location as a function of £ and MEdd 

For a particular value of £, MEdd and 7i„, the shock location (measured from the black hole in units of r g ) can be 
calculated using Eq. (10). As u s h and M s h is a function of £, MEdd, and ji n , r s h will also change with the change of 
any of these accretion parameters. In Fig.l, we show the variation of r s with £ and MEdd for a fixed j in — |. While 
£ and MEdd are plotted along X and Y axes respectively, r s h is plotted along Z axis in units of r g . It is interesting to 
note that different curves start and terminate at different points, which indicates that shock formation is not a generic 
phenomena, i.e., shock does not form for any value oi£ and MEdd', rather, a specific region of parameter space spanned 
by jf , MEdd, 7m| allows shock formation. Also wc find that depending on the value of £ , both sub-Eddington as well 
as super-Eddington accretion allows shock formation for a fixed value of 7$ n . It is evident from the figure that while 
shock location non- linearly anti-correlates with £ (for a given value of MEdd and 7m), it correlates with MEdd (for 
a fixed £ and 7»„). Also we see that the maximum value of energy £mo,x f° r which the shock forms decreases with 
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Fig. 2: Combined solution topology for transonic accretion-wind system for fixed values (shown in the figure) of accretion 
parameters S, MEdd, Jin and for four outflow polytropic indices y — 1.3, 1.275, 1.25 and 1.225 respectively (from left to 
right). Solid curve marked with AB represents the pre-shock transonic accretion while four dotted curves represent the accretion- 
powered outflow branches. Solid vertical line BC (marked by r s h) with double arrow stands for the shock transition and solid 
line marked by CD stands for the 'self-wind' branch. P a is the location of the inflow sonic point. See text for details. 



an increase of MEdd and also the maximum as well as minimum value of MEdd for which the shock forms, decreases 
and increases respectively with an increase in £ for a fixed -f in . Families of curves for other values of j in can also be 
obtained in the same manner . 



3.2. Combined integral curves of motion 

Fig. 2 shows a typical solution which combines the accretion and the outflow. While Mach number is plotted along 
the Y axis, the distance (in units of r g ) from the event horizon of the accreting black hole is plotted along the X 
axis in logarithmic scale. The accretion parameters used are £ = 1.0008, M-Edd =0.5 and 7j„ = | corresponding to 
ultra-relativistic inflow. The solid curve AB represents the pre-shock region of the inflow and the solid vertical line 
BC with double arrow at r s h represents the shock transition. Shock location (3.8 r g ) is obtained using the Eq.(10) for 
a particular set of inflow parameters mentioned above. Four dotted curves show the four different outflow branches 
corresponding to the four different adiabatic indices j of the outflow. From left to right, the values of 7„ are 1.3, 1.275, 
1.25 and 1.225 respectively with respective mass-outflow rates as 0.0002, 0.00006, 0.000011 and 0.00000122 respectively 
which indicates that for a given value of £, MEdd and j in , R in correlates with j , which is explicitly shown in Fig 
6. It is evident from the figure that the outflow moves along the solution curves in a completely different way to 
that of the 'self- wind' solution of the inflow (solid line marked by CD). Also, the sonic points for all the outflowing 
branches are different to those of the accretion 'self-wind' system which is designated as P s . While P s = 470.8 r g , the 
sonic points of the outflowing branches corresponding to j = 1.3, 1.275, 1.25 and 1.225 are 574.96, 669.7, 783.34 and 
922.3 r g respectively, which indicates that the outflow sonic point increases with a decrease in the adiabatic index of 
the outflow and thus the wind starts with a very low bulk velocity which is why the mass-loss rate decreases. It is 
also observed that the sonic point of the accretion-'self-wind' system is, in general, located closer to the event horizon 
compared to the outflow sonic point for all values of £, MEdd, Jin and j Q . 
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Fig. 3: Variation of Rm with conserved total specific energy (which includes the rest mass energy) £ (dotted curve in the figure 
marked with £) and post shock temperature (solid curve marked with T^ h ) of the flow for a fixed value of accretion rate and 
polytropic indices of the flow having values shown in the figure. Both £ and Tf h are plotted along X axis with the following 
scaling: 

(where Tn = 10 11 °K), to fit in the same graph. 



1 sh 



1 + 0.0004 




Fig. 4: Variation of Rm with accretion rate (dotted line marked MEdd) and corresponding post-shock flow temperature (solid 
line in the figure marked by T^ h ) for fixed values of £ , ji n , and 7 shown in the figure. T^ h is scaled as T^ h — > 2.2 • If 

is evident from the figure that low-luminosity objects produce more mass-loss, see text for details. 



3.3. Dependence of R m on various flow parameters 
3.3.1. Variation of R m with £ 

In Fig. 3, we have plotted the variation of with conserved inflow specific energy (which includes the rest mass 
energy of matter) £ (dotted line marked with £) for a fixed value of MEdd (= 0.5), 7i„ (= |) and 7 Q (= 1.3). We 
observe that the mass-outflow rate non-linearly correlates with £. The explanation may be as follows: 
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Fig. 5a: For fixed values of energy and accretion rate shown in the figure, variation of Rm with polytropic index of the inflow 
(7»n) is presented for five different values of A 7 = 7i„ — 7 , where 7 is the polytropic index of the outflow. It is observed that 
ultra-relativistic flows produce less mass-loss compared to pure-nonrelativistic flows, see text for details. 
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Fig. 5b: Variation of Rm with post-shock flow temperature T^ h corresponding to values of 7^ 
is observed to be thermally driven, see text for detail. 



shown in the figure 5(a). Outflow 



As £ increases, r s h decreases (see Fig. 1, §3.1) and the post-shock bulk velocity of the flow u^ h as well as the post- 
shock density p~^ h increases. The outflow rate, which is the product of three quantities r s h, p^ h and u^ h (see Eq. (15)), 
increases in general due to the combined 'tug of war' of these three quantities. Moreover, the closer the shock is to the 
black hole, the greater the amount of gravitational potential will be available to be put onto the relativistic protons 
to provide stronger outward pressure and the closer the shock forms to the black hole, the higher is the post-shock 
temperature (the effective characteristic outflow temperature) and the higher is the amount of outflow (as the wind is 
observed to be strongly thermally driven, see discussion below). Thus the mass-outflow rate increases with £ because 
for a particular set of fixed values of MEdd, Rm is proportional to M out which increases with £. Also in the same 
figure, we show the variation of R m with post-shock temperature T^~ h (solid line in the figure marked with T^ h ). The 

value of T^ h corresponds to the values of £ shown in the same figure and is scaled as T^ h — ► 1 + 0.0004 x (7^7) 

(where Tn = 10 11 °K) to fit in the same graph. We see that post-shock temperature correlates with the energy of the 
flow, and for a fixed accretion rate and adiabatic indices of the inflow and outflow, R m also correlates with post-shock 
temperature, which indicates that outflow is thermally driven as well. If we draw a family of curves for R m vs. £ 
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Fig. 6: Variation of Rm with polytropic index of the outflow j for a fixed values of £ , M B dd and ~f in as shown in the figure. 

(alternatively, R m vs. T^ h ) for different values of MEdd (figure not presented in this paper) ranging from sub-Eddington 
to super-Eddington accretion rate, we observe that the maximum value of R m obtained for a fixed MEdd (while £ 
is being varied) increases with a decrease of MEdd and also the maximum value of T^ h follows the same trend as 
mentioned above. So a 'high energy-low luminosity' combination for polytropic accretion maximizes the post-shock 
flow temperature and gives rise to the highest amount of outflow. 

3.3.2. Variation of R m with Mem 

In Fig. 4 we show the variation of R m as a function of accretion rate (dotted curve marked with MEdd) scaled in units 
of Eddington rate for a fixed value of £ (=1.00045) ji n (= |) and -f Q (=1.3). It is observed that R m anticorrelates 
with MEdd which indicates that low-luminosity flow produces more mass-loss. Here R m correlates with post shock 
temperature T^ h (solid curve marked with T^ h ). T^ h shown here is scaled as T^ h — ► 2.2 x (jr^j to fit in the same 

graph. The values of T^ h correspond to the varying MEdd shown in this graph. It is also obvious from this figure that 
outflow is thermally driven because R m correlates with T^~ h . It can be shown that while for a fixed value of £ and 
lin-, r sh,pth an d u th correlates with MEdd, shock Mack number, compression ratio and post-shock flow temperature 
anti-correlates with the accretion rate. However, from the figure we see that unlike £ , variation of R m is somewhat 
insensitive to the variation of MEdd for fixed £ and ji n . 

3.3.3. Variation of R m with adiabatic indices 

In previous cases, the polytropic index of the accreting matter was always kept fixed at the value |. To have a 
better insight of the behaviour of the outflow, we plot Rm as a function of 7j„ (fig 5a) for a fixed value of £ (= 1.0002) 
and MEdd (= 0.2). The upper range of 7j„ shown here is the range for which shock forms for the specified value of £ 
and MEdd- We have chosen the value of j Q in such a way so that the condition j a < j in is always satisfied. Defining 
A 7 to be A 7 = ji n — 7 , we study the variation of R m with -fi n for five values of A 7 = 0.01, 0.02, 0.03, 0.04 and 0.05 
(from top to bottom respectively). We observe that R m correlates with -f in , which is expected because the specific 
enthalpy of the flow increases with j in to produce a higher post-shock temperature for higher value of j in (see also 
Fig. 5(b)). We observe from our calculation that p~^ h and u s h correlates while r s h anticorrelates with 7j„. So increment 
of "fin satisfies all possible conditions to have a high value of R m . In fig. 5b, we show the variation of R m with Tf h 
(corresponding to the values of 7j„ shown in fig. 5a) scaled in units of Tn = 10 10 °K to show that here also the flow is 
thermally driven. Also, we observe that M a h as well as the shock compression ratio R CO mp anti correlates with 7j„ so 
that 'weak-shock' solutions are preferred to obtain a high value of mass-loss for this case. So we conclude that as the 
accretion approaches from its ultra-relativistic nature to its non-relativistic regime, the mass-loss rate increases. 
In Fig. 6, we show the variation of R m with 7 Q for a fixed value of £ = 1.0005, MEdd = 0.2 and ji n = |. The general 
conclusion is that R m correlates with 7 . This is because as -f Q increases, shock location and post-shock density of 



i i | i i i | i i i | i i i | r 

5=1.0005 




Das T. K.: Accretion Powered Spherical Wind in General Relativity 



13 



matter does not change (as "f does not have any role in shock formation or in determining the R CO mp) but the sonic 
point of the outflow is pushed inward, hence the velocity with which outflow leaves the shock surface goes up, resulting 
the increment in R m . 

4. Conclusion 

In this paper, we could successfully construct a self-consistent spherically-symmetric, polytropic, transonic, non- 
magnetized inflow-outflow system by simultaneously solving the set of hydrodynamic equations governing the ac- 
cretion and wind around a Schwarzschild black hole using full general relativistic framework. Introducing a steady, 
standing, hadronic-pressure supported spherical shock (formation of which was first proposed by KE86 and PK83) 
surface around the black hole as the effective physical atmosphere which may be responsible for generation of accretion- 
powered spherical wind, we calculate the mass-outflow rate R m in terms of only three accretion parameters (conserved 
energy of the flow £ , which includes the rest mass energy of matter, accretion rate Msdd scaled in units of Eddington 
rate and polytropic index of the flow j in ) and only one outflow parameter (the polytropic index of the outflow, j a ). 
Not only do we provide a sufficiently plausible estimation of Rm, we could also successfully study the dependence and 
variation of this rate on various physical parameters governing the flow. 

At this point, it is worth mentioning that the hot and dense shock surface around black holes, which is proposed here 
as the effective physical barrier around compact objects regarding the mass outflow, may be generated due to other 
physical effects as well for spherical accretion (Chang & Osttriker 1985, Meszaros & Ostriker 1983, Babul, Ostriker & 
Meszaros 1989, Park 1990, 1990a). Another very important approach launched recently was to construct such an 'ef- 
fective barrier' for non-spherical disc accretion to introduce the concept of CENtrifugal pressure supported BOundary 
Layers (CENBOL). Treating the CENBOL as the effective atmosphere of the rotating flows around compact objects 
(which forms as a result of standing Rankine-Hugoniot shock or due to the maximization of polytropic pressure of 
accreting material in absence of shock), detailed computation of the mass-outflow rate from the advective accretion 
disks has been done, and dependence of this rate on various accretion and shock parameters has been quantitatively 
studied by constructing a self-consistent disk-outflow system (Das 1998, Das & Chakrabarti 1999). 
The basic conclusions of this paper are the following: 

1. Shock formation is not a generic phenomena, i.e., not all solutions contain shock, rather a specific region of 
parameter space spanned by £, MEdd and ji n allows shock formation. Also we found that for given values of 
£,MEdd and j in , while the value of shock location (in units of r g ) correlates with MEdd, it anti- correlates with 
both £ and ji n . 

2. The shock surface can serve as the 'effective' physical barrier around the black hole regarding generation of mass 
loss via transonic spherical wind. The fraction of accreting material being blown as wind (which is denoted as R m ) 
could be computed in terms of three accretion parameters and one outflow parameter. 

3. While R m correlates with £,%„ and it anti- correlates with MEdd, which indicates that low-luminosity objects 
produce more mass-loss, though outflow could be generated for both sub-Eddington as well as super-Eddington 
accretion. 

4. If a shock forms, then whatever the initial flow conditions and whatever the nature of dependence of R m on any of 
the accretion/ shock/ outflow parameters, R m always correlates with post-shock flow temperature, which indicates 
that outflow is strongly thermally driven; hotter flow always produces more winds. 

Our calculations in this paper, being simply founded, do not explicitly include various radiation losses and cooling 
processes, combined effects of which may reduce the post-shock proton temperature (which means the reduction of 
outflow temperature), in reality could be lower than what we have obtained here and the amount of outflow would be 
less than what is obtained in our calculation. This deviation will be more important for systems with high accretion 
rates. Nevertheless, cases of low accretion rates discussed here would not be affected that much and our preliminary 
investigation shows that even if we incorporate various losses, the overall profile of the various curves showing the 
dependence of R m on different inflow parameters would be exactly the same, only the numerical value of R m in some 
cases (especially for high accretion) might decrease. However, as we have shown that R m anti-correlate with accretion 
rate ( see §3.3.2, Fig. 4) and also it has been shown that the mass-loss rate does not have that strong a dependence 
on Msdd, we believe that such changes will not bring any significant change in our overall conclusions. 
Our present work, as we believe, may have some important consequences regarding modelling of the astrophysical 
outflows emanating from galactic and extra-galactic sources powered by isolated accreting compact objects. Galactic 
and extra-galactic sources of jets and outflows are now widely believed to be fed by accreting black holes sitting at 
the dynamical centre of these sources. In the absence of any binary companion, spherically symmetric Bondi (Bondi 
1952) type accretion may occur onto isolated central black hole if accreting matter has a negligible amount of intrinsic 
angular momentum. This may happen when the central super-massive black hole at the galactic centre is surrounded 
by a dense stellar cloud in such a way that the vector sum of the angular momentum of tidally disrupted matter (from 
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a number of stars with trajectories approaching sufficiently close to the hole) almost vanishes. On the other hand, 
unlike the ordinary stellar bodies, black holes do not have their own 'physical' atmosphere and outflows in these cases 
have to be generated from the accreting matter only. Hence we believe that it is necessary to study the accretion and 
outflow (from various astrophysical sources powered by accreting compact objects) in the same framework instead of 
treating the outflow separately from accretion phenomenon. At the same time, as the fundamental criterion for con- 
structing any self-consistent physical model demands the minimization of the number of inputs to the model, we may 
conclude that modelling the astrophysical outflow needs a concrete formulation where the outflow can be described 
in terms of minimum number of physical parameters governing the inflow. The success of our work presented in this 
paper, as we believe, is precisely this. We could rigorously compute the mass outflow rate R m in terms of only three 
inflow parameters and only one outflow parameter. To the best of our knowledge, there is no such model present in the 
literature which rigorously studies the accretion-powered spherical outflow in a general relativistic framework using 
only the four parameters mentioned above. 

Although in this work we have performed our calculation for a I 0M© Schwarzschild black hole, general flow charac- 
teristics will be unchanged for black hole of any mass except the fact that the region of parameter space responsible 
for shock formation will be shifted and the value of R m will explicitly depend on the mass of the black hole. In our 
next work, we will show the direct dependence of R m on the mass of the black hole and will apply our model to some 
specific astrophysical sources of outflows to roughly estimate the amount of outflow in M Q /j/r _1 . 
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